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Abstract 

Consider a directed or an undirected graph with integral edge weights from the set [— W, W], 
that does not contain negative weight cycles. In this paper, we introduce a general framework 
for solving problems on such graphs using matrix multiplication. The framework is based on 
the usage of Baur-Strassen's theorem and of Strojohann's determinant algorithm. It allows us 
to give new and simple solutions to the following problems: 

Finding Shortest Cycles - We give a simple O(Wn 0J ) time algorithm for finding shortest cy- 
cles in undirected and directed graphs. For directed graphs (and undirected graphs with non- 
negative weights) this matches the time bounds obtained in 2011 by Roditty and Vassilevska- 
Williams. On the other hand, no algorithm working in 0{Wn LJ ) time was previously known for 
undirected graphs with negative weights. Furthermore our algorithm for a given directed or 
undirected graph detects whether it contains a negative weight cycle within the same running 
time. 

Computing Diameter and Radius - We give a simple 0(Wn u ) time algorithm for comput- 
ing a diameter and radius of an undirected or directed graphs. To the best of our knowledge 
no algorithm with this running time was known for undirected graphs with negative weights. 

Finding Minimum Weight Perfect Matchings - We present an 0(Wn u ) time algorithm 
for finding minimum weight perfect matchings in undirected graphs. This resolves an open 
problem posted by Sankowski in 2006, who presented such an algorithm but only in the case 
of bipartite graphs. 

We believe that the presented framework can find applications for solving larger spectra of re- 
lated problems. As an illustrative example we apply it to the problem of computing a set of 
vertices that lie on cycles of length at most t, for some t. We give a simple 0{Wn") time algo- 
rithm for this problem that improves over the 0(Wn^t) time algorithm given by Yustcr in 2011. 

In order to solve minimum weight perfect matching problem we develop a novel combinatorial 
interpretation of the dual solution which sheds new light on this problem. Such a combinatorial 
interpretation was not know previously, and is of independent interest. 
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1 Introduction 



The application of matrix multiplication to graph problems has been actively studied in recent 
years. The special case of unweighted graphs is well understood. For example, 0{n w ) time al- 
gorithms for finding shortest cycles [16] have been known for 35 yeara^J But similar results for 
weighted graphs were obtained only last year, by Roditty and Vassilevska- Williams [25J. Their 
algorithm works in OiWn^) time, where W is the largest magnitude of an edge weight. Two sim- 
ilar problems on weighted graphs where there has been considerable effort, but the full answer has 
not been achieved, are diameter and perfect weighted matching. This paper introduces a general 
framework that gives simpk^] solutions to all three of these problems and others. In the following 
unless otherwise stated, we work with graphs that contain edges with possibly negative weights but 
no negative cycles. We obtain the following results. 

Finding Shortest Cycles We give a simple 0{Wn UJ ) time algorithms for finding shortest cycles 
in undirected and directed graphs. In the case of directed graphs the algorithm reduces the problem 
to one determinant computation for a polynomial matrix. On the other hand, the undirected case 
requires handling short 2-edge cycles in a proper way. The idea used here is an extension of the 
algorithm by Sankowski contained in (27], that allowed to test whether a graph contains negative 
weight cycle. For directed graphs (and undirected graphs with non- negative weights), our bounds 
match the ones obtained in 2011 by Roditty and Vassilevska- Williams [25], whereas for undirected 
graphs with negative weights no 0{Wn u ) time algorithm was previously know for this problem. 
For all related results see Table 1 and Table 2. For undirected graphs with negative weights the 
problem reduces to n computations of minimum weight perfect matchings [6]. 

Furthermore our algorithm for a given directed or undirected graph detects whether it contains 
a negative weight cycle within the same running time. 

Computing Diameter and Radius We present a simple 0{Wn u> ) time algorithm for comput- 
ing a diameter and radius of undirected and directed graphs. This algorithm combines determinant 
computations with binary search. Since computing all pairs shortest paths suffices to find both 
diameter and radius 0(Wn UJ ) time algorithm follows from [31] in case of undirected graphs with 
non-negative weights. Moreover by generalizing the ideas of [31] used for non- negative weights by 
applying random sampling one can obtain the same running time for directed graphs without neg- 
ative weight cycles. However, to the best of our knowledge, all previous solutions to this problem 
in undirected graphs with negative weights reduced the problem to n computations of minimum 
weight perfect matchings. For other related results see Table 3 and Table 2. 

Finding Minimum Weight Perfect Matchings We present an 0(Wn u ) time algorithm for 
finding minimum weight perfect matchings in undirected graphs. This resolves an open problem 
posted by Sankowski in 2006 [28J , who presented such an algorithm but only in the case of bipartite 
graphs. Some advance on this problem has been recently given by Huang and Kavitha [15] . who 
have shown an 0{Wn u ) time algorithm for the maximum weight matching problem. However, the 

1 The O notation ignores factors of log n and logW. O(n^) is the time needed for a straight-line program to 
multiply two n x n matrices; cj is called matrix multiplication exponent. 

2 An objective sense in which our algorithms are simple is their use of algebra: The power of our algebraic 
algorithms comes from black-box routines, and the algorithms themselves use only elementary algebraic ideas. 
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weighted perfect matching problem is more involved and no reduction similar to the one presented 
in [15] is known to work. Previously, similar reduction was given in [18] for maximum weight 
bipartite matching. Nevertheless, to solve minimum weight perfect matching problem even in bi- 
partite graphs one is bound to use more structured techniques. Actually, we need to develop such 
a technique for general graphs. We give a novel combinatorial interpretation of the dual problem. 
Such an interpretation for general matching problem was not know previously and is of significant 
independent interest. For the summary of different algorithms see Table 2. 



Complexity 


Author 


0(nm + n 2 log n) 
directed 


Johnson (1977) 117| 


0(n") nonnegative un- 
weighted 


Itai and Rodeh (1977) 


0(VK u b81 n 2 bYb ) directed 


Zwick (2000) [37| 


0(nm + n 2 log log n) directed 


Pettie (2004) 24 


0(n 3 log 15 logn/ log 2 n) 
directed 


Chan (2007) g] 


0(Wn hJ ) directed and 
nonnegative undirected 


Roditty and Vassilevska- 
Williams (2011) [25] 


0(Wn") 


this paper 



Table 1: The complexity results for the shortest cycle 
problem. 



Complexity 


Author 


0(nm + n 2 logn) directed 


Johnson (1977) [171 


0(n min(m log n, n 2 )) undi- 
rected 


Gabow (1983) [8] 


0(Wn UJ ) nonnegative undi- 
rected 


Shoshan and Zwick 
(1999) [31] 


O(VK b81 n 2 bYb ) directed 


Zwick (2000) 137J 


0(n :i log Li logn/ log 2 n) 
directed 


Chan (2007) 0] 


0(Wn UJ ) directed 


folklore+|311+[37| 


0(Wn u ) 


this paper 



Table 3: The complexity results for the diameter and 
radius problem. Note that the first five results solve 
the more general All Pairs Shortest Paths problem. 



Complexity 


Author 


0(n 2 m) 


Edmonds (1965) [5] 


0(n 3 ) 


Lawlcr (1973) [20J and 
Gabow (1974) [7J 


0(nm log n) 


Galil, Micali and 
Gabow (1982) Q3] 


0(n(m log log log m n+nlogn)) 


Gabow, Galil and 
Spencer (1984) [11] 


Ofnl m log W) 


Gabow (1985) [9] 


0(n(m + n log n)) 


Gabow (1990) [10] 


0(^/na(m, n) logn mlog(nW)) 


Gabow and Tarjan 
(1991) [12] 


0(Wn") 


this paper 



Table 2: The complexity results for the minimum 
weight perfect matching problem. Note that first, 
third, fourth and fifth results solve the more general 
All Pairs Shortest Paths problem. 



Complexity 


Author 


0(nm) 


Suurballc and Tarjan (1984) 




[33] 


0(tWn") undirected 


Yuster (2011) [34] 


0(t 4 - u Wn u ) directed 


Yuster (2011) [34] 


O(W) 


this paper 



Table 4: The complexity results for the problem of 
finding vertices that lie on cycles of length at most t. 



1.1 The Framework 

Our framework is based on two seminal results. First of all we use Storjohann's algorithm |32j that 
computes the determinant of a degree d polynomial matrix in 0{dn UJ ) time. All the above graph 
problems can be encoded as a determinant problem on a polynomial matrix. However the determi- 
nant is just one number and does not provide enough information to decode the whole solution. Here 
we use the Baur-Strassen Theorem [2], which shows how to compute all partial derivatives of a func- 
tion in the same asymptotic time as computing the function itself. For a simple constructive proof 
please check [22J . This theorem allows us to magnify the output of the algorithm from 1 number to 
n 2 numbers. The algorithms obtained in this way are very simple and work in three phases: compute 
the determinant of an appropriately defined matrix; apply Baur-Strassen to the result; decode the 
output. Even for minimum weight matching our algorithm is simple, and computes the dual solution 
in just a few lines of pseudocode. Based on these examples, we believe this framework will find appli- 
cation to a wider spectrum of problems. The full version of the paper will contain a more extensive 
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review of possible applications. Here we give one more illustration: a simple 0(Wn u ) time algorithm 
that finds every vertex that lies on a cycle of length < t, for a given arbitrary t (for a directed or undi- 
rected weighted graph, with negative edges allowed but no negative cycles). This improves the algo- 
rithms proposed by Yuster in 2011 [M] (see Table 4; those algorithms do not allow negative edges). 

The paper is organized as follows. In the next two sections we introduce needed tools and give 
the main definitions. In Section [4] we introduce our framework and give algorithms for the shortest 
cycle problem in directed graphs. This motivates the introduction of the framework in Section [5| 
Section [6] contains the algorithm for minimum weight perfect matching problem. The next section 
presents the application of our framework for computing the diameter of a graph. Section [8] intro- 
duces the ideas needed to compute shortest cycles in undirected graphs with positive weights. In 
Section [9] we join the ideas from all previous section to solve the shortest cycle problem and the 
diameter problem in undirected graphs with negative weights. Finally, in Section 10 we apply our 
framework to find the set of vertices that lie on a cycle of lenght at most t. 



2 Preliminaries 

Our approach is based on three main ingredients. 



Linear Algebra Algorithms Storjohann [32] has made an important addition to the set of prob- 
lems solvable in 0(n u ) arithmetic operations: the determinant and the rational system solution for 
polynomial matrices. 

Theorem 1 (Storjohann '03). Let K be an arbitrary field, A G K[y] nxn a polynomial matrix of 
degree d, and b £ K[y] nxl a polynomial vector of the same degree. Then 

• rational system solution A~ l b (Algorithm 5 132]), 

• determinant det(A) (Algorithm 12 [32]), 

can be computed in 0{dn u ) operations in K, with high probability. 

For the next section, note that both algorithms of Storjohann can be written as straight-line 
programs. The randomization does not pose a problem as it is just used at the very beginning 
to generate a polynomial of degree d that does not divide det(^4). After this the algorithms are 
deterministic. Our applications will work over a finite field (see Section [2]) so there is no risk 
of manipulating huge integers. Finally usage of the FFT to perform multiplication of degree nd 
polynomials does not pose a problem. 



Baur-Strassen Theorem Another astonishing result is the Baur-Strassen Theorem from 1983 [2, 
122] . It was used to show that matrix multiplication is no harder than determinant computation 
when considering algorithms that can be written as straight-line programs. Such a statement is 
unexpected when you realize that matrix multiplication returns n 2 numbers and determinant com- 
putation just one. However it is possible to increase the number of outputs by modifying the 
algorithm appropriately. Let T(/i, ...,/&) denote the time needed to compute functions /i, . . . , fk, 
all at the same given point. 

Theorem 2 (Baur-Strassen '83). For straight-line 

T(/,f^,...,^)<5T(/). 
OX i ox n 
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Thus we can compute all partial derivatives of a function / in the same asymptotic time as is 
used to compute /. [22] shows a RAM routine to compute all n partials can be constructed in time 
0(T(/)) as well. 



Schwartz-Zippel Lemma Our approach is to encode a graph problem in a symbolic matrix 
whose determinant is (symbolically) non-zero if and only if the problem has a solution. The 
Schwartz-Zippel Lemma [36, 30] provides the non-zero test: 



Corollary 3. For any prime p, if a (non-zero) multivariate polynomial of degree d over Z p is 
evaluated at a random point, the probability of false zero is < d/p. 

We will choose primes p of size B(n c ) for some constant c. In a RAM machine with word size 
G(logra), arithmetic modulo p can be realized in constant time. Moreover we would like to note 
that multivariate determinant expressions have been used several times in previous work, including 
Kirchhoff 's Matrix- Tree Theorem, results of Tutte and Edmonds on perfect matchings, and recent 
result of Bjorklund [3] for undirected Hamiltonicity. 

3 Definitions: Graphs with Integral Weights 

A weighted n-vertex graph G is a tuple G = (V, E,w,W), where the vertex set is given by 
V = {1, ... , n}, E C V x V denotes the edge set, and the function w : E — >• [— W, W] ascribes 
weights to the edges. In this paper we consider both undirected and directed graphs. Hence, E 
might denote either a set of unordered pairs for undirected graphs or a set of ordered pairs for 
directed graphs. We define a weight of the edge set F C E to be e(F) = ^2 e£ p w(e). 

Consider a path p = v\, V2, ■ ■ ■ , of length k. The weight of this path is simply equal to the 
weight of the edge set of p and denoted by by w(p). The distance from v to u in G, denoted by 
distc(v,u), is equal to the minimum weight of the paths starting at v and ending in w. A path 
from v to u with minimum weight is called a shortest path. If for a path p we have v\ = then 
the path is called a cycle. 

In the shortest cycle problem we are given a weighted graph G and we need to compute the 
shortest (i.e., minimum weight) cycle in G. In this problem we assume that the graph contains no 
negative weight cycles. There is a subtle difference between the directed and undirected versions 
of the problem. The standard approach to reduce the undirected problem to directed problem by 
bidirecting the edges does not work. The resulting graph can contain cycles that pass though the 
same edge twice in both directions. Such cycles do not exist in the undirected graph. Moreover, 
when the undirected graph contains negative edges the resulting directed graph contains negative 
weight cycles, even if the undirected graph did not. 

In the diameter problem we are given a weighted graph G (directed or undirected) and are asked 
to find the a pair of vertices v,u £ V that maximizes distG(v,u). Note that here the reduction 
of undirected problem to directed one by bidirecting the edges works on the in the case when the 
edges are positive. 

A matching M in G is a set of edges such that each vertex is incident at most one edge from 
M . In a perfect matching each vertex is incident to exactly one matched edge. A minimum weight 
perfect matching is a perfect matching M in a weighted undirected graph that minimizes the total 
weight w{M) = ^egM w ( e )- Many other notions of "optimum weighted matching" reduce to min- 
imum weight perfect matching: A maximum weight perfect matching is equivalent to a minimum 
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weight perfect matching for weights w'(e) := — w(e). A minimum weight cardinality k matching 
(i.e., exactly k edges are to be matched) is a minimum weight perfect matching on the graph with 
n — 2k articial vertices, each joined to every original vertex by a zero-weight edge. A maximum 
weight matching (i.e., we want to maximize the total weight of the matched edges) is a maximum 
weight perfect matching on the graph that has one artificial vertex if n is odd, plus new zero-weight 
edges that make the graph complete. There is a reduction going in the other direction: we set edge 
weights to w'(e) := nW + w(e). However we are interested in time bounds that are linear in W, 
so this reduction is not of interest. 



4 Shortest Cycles in Directed Graphs 

Let K be an arbitrary field. A symbolic polynomial p[y] is a multivariate polynomial over a set 
of variables y U X over K. We denote the set of symbolic polynomials by K[y] = K[y U X]. We 
write as well K to denote the set of multivariate polynomials i^[A]. A symbolic matrix polynomial 
A[y] £ K[y] nxn is an n x n matrix whose entries are symbolic polynomials from K[y]. 

We shall use a straight-line program that evaluates det(yl[y]) {A[y\ £ K[y] nxn ) using Storjo- 
hann's algorithm. Here the goal is to evaluate the determinant to a polynomial in one variable 
y. This program is easily constructed: Start with the original straight-line program that evaluates 
det(^4) (A £ K[y] nxn ) using Storjohann's algorithm. Prepend assignment statements of the form 
a-ijk ^— (iijk, where is the variable in the original program for the coefficient of y k in the en- 
try Aij, and the corresponding coefficient in A[y] is G K. In our applications each of these 
new assignment statements uses O(l) time, so the extra time can be ignored. Also note that all 
arithmetic in our straight-line programs is done in K = Z p for a chosen prime p. 

We say that a : X — > K is an evaluation function. We define define to be a one variable 

polynomial in y with all variables i£l substituted by a(x). 

Let us define a symbolic adjacency matrix of the weighted directed graph = (V, E, w, W) to 
be the symbolic matrix polynomial A^) such that 



' J \ otherwise, 



where Xi j are unique variables corresponding to the edges of Hence, X is the set of all variables 
Xij. For a multi- variable polynomial q, let us denote by: 

• deg*(q) - the degree of the smallest degree term of y in q, 

• term^(g) - the coefficient of y d in q, 

• term*(c/) - term^(^) for d = deg*(q). 

If q is the polynomial, deg*(q) := oo. In the following we assume that K := Z p , i.e., we work 
over a finite field of order p for some prime number p. 

We say that a nonempty set of disjoint cycles C in ^ is a cycle packing if every vertex of 
belongs to at most one cycle in C. Observe that a minimum weight cycle packing is either a shortest 
cycle or a collection of shortest cycles all of weight 0. 

Lemma 4. Let ~di be a directed weighted graph. The minimum weight of a cycle packing in 
equals deg*(det(^L(^) +/) — !). This weight is also the weight of a shortest cycle if G has no 
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negative cycle. Finally for any G and any p > 2, all non-zero terms in det(A(G) +1) are non-zero 
over a finite field Z p . 

Proof. By definition 

n 

det (A0) +/) = E ff (P) II { A $) +l) k . (1) 

pe r„ fe=i ,Pfc 

where r„ is the set of n-element permutations, and cr(p) is the sign of the permutation p. A permu- 
tation p defines a set of directed edges C p = {(i,Pi) : 1 < i < n}. This edge set corresponds to a set of 
cycles given by the cycles of p. The edge set C p includes self- loops for all i such that i = Pi- Note that 
p corresponds to a non-zero term in the determinant if and only if C p contains only edges from E or 
self-loops. Hence, after throwing away self-loops p can be identified with a cycle packing in G, or 0. 

Let us now compute the degree of y in the term of the determinant given by p. This term is 
obtained by multiplying the elements of the matrix A(G) + I corresponding to the edges of C p . The 
degree of the term equals the sum of the degrees in its individual elements. These degrees encode 
the weights of the edges or are zero for self-loops. So the degree of the term is the total weight of 
the cycles in C p , where self- loops have weight zero. 

The term corresponding to a set of self-loops is equal to 1. Hence, det(A(G) + 1) — 1 contains 
only terms that correspond to nonempty sets of cycles. So we have proved the first assertion, i.e., 
the smallest degree of y is the smallest weight of a packings of cycles. 

For the third assertion note that each monomial in ([!]) has coefficient ±1 as each of them 
contains different variables. Hence each non-zero term is also non-zero over Z p for any p > 2. 

Finally for the second assertion assume G has no negative cycle. Let C be the shortest cycle 
in It corresponds to a cycle packing of weight w(C). (The packing has a self-loop for every 
v C.) In fact this is the minimum weight cycle packing. In proof take any cycle packing. If it 
contains more than one non-loop cycle discard all but one of them to get a packing of no greater 
weight (since every cycle is nonnegative). By definition this weight is > w(C). □ 

Computing the Weight of the Shortest Cycle Next we apply the ideas of the previous sec- 
tion to compute the weight of a shortest cycle. Since A{G) may have entries with negative powers 
of y, we cannot use Storjohann's result to compute its determinant. However we can multiply it 
by y w to make exponents nonnegative and use this identity: 

de g ;(det(i(^) + I) - 1) = deg*(det((i(^) + I)y w ) - y nW ) - nW. 
Combining this idea with the Schwartz-Zippel lemma we obtain the following algorithm. 

Algorithm 1 Computes the weight of the shortest cycle in a directed graph G. 
1: Generate a random substitution a : X — > Z p for a prime p of order 0(n 2 ). 
2: Compute 5 = det((A(G*)\ a + I)y w ) — y nW using Storjohann's theorem. 
3: Return deg*(5) — nW . 



This algorithm implies. 
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Theorem 5. Let G = (V, E,w,W) be a weighted directed graph without negative weight cycles. 
The weight of the shortest cycle in ~G can be computed in 0{Wn^) time, with high probability. 

Note the algorithm also detects the existence of a negative cycle- it corresponds to 5 < 0. 

Finding a Shortest Cycle The above algorithms does not seem to give the cycle by itself. There 
is, however, a straightforward way of doing it using the Baur-Strassen theorem. We say that edge 
e is allowed if and only if it belongs to some shortest cycle C in G. 

Lemma 6. The edge (u, v) 6 E is allowed if and only if: term* (det(A((f )+/) — 1) is non-zero. 
Moreover, for p > 2, it is non-zero over a finite field Z p . 

Proof. The above follows directly by the proof of Lemma [4} The partial derivative is non-zero if 
and only if the variable x U)V exists in some smallest degree term, and so the corresponding edge 
(u, v) has to lie on some shortest cycle. 

□ 

Rewrite the expression of the lemma to eliminate negative powers of y: 



term*,(det(I(<^) + /)-!) = ter m ;(det((10) + /) 



This implies the following algorithm and theorem. 



V ) - V 



Algorithm 2 Computes a shortest cycle in a directed graph C^. 



1: Let d x f be the routine given by the Baur-Strassen Theorem to compute the matrix of partial 

is the degree deg*(<5) computed 



derivatives q® term^ 
in Algorithm 1. 

Generate a random substitution a : X — > Z p for a prime p of order 0(n 2 ). 

Compute the matrix 5 = d x f\ a . 

Take any edge (u, v) such that 8 U , V ^ 0. 

Compute the shortest path p V:U from v to u in using [27] [35] . 
Return the cycle formed by (n, v) and p VjU . 



Note that in Step 1 we are applying the Baur-Strassen Theorem to the straight-line program 



forterm^ det{(A(G) + I)y ) — y ■ The latter is constructed using Storjohann's Theorem and 
the modification described at the start of this section. 

The shortest path computation using algorithms from [27}I35| takes O^Wn^) time and succeeds 
with high probability, so we obtain: 

Theorem 7. Let = (V, E, w, W) be a weighted directed graph without negative weight cycles. 
The shortest cycle in can be found in OiWn w ) time, with high probability. 
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5 Performing Symbolic Computations 



This section gives a formal description of the algebraic operations used in our algorithms. Let A[y] 
be an n x n symbolic matrix and let a : X — > Z p be an arbitrary evaluation function. Moreover for 
a prime p let / : Z p [y] nxn — > Z p be a symbolic matrix function. 

Theorem 8. If the result of f has degree bounded by some polynomial dj(n) and f(A[y]\ a ) is com- 
putable in tf(n) time, then there exists an algorithm that in 0(tf(n)) time checks whether f(A[y]) 
is symbolically non-zero over Z p with error probability < df(n)/p. 

Proof. The algorithm works as follows: 

1. uniformly at random choose a substitution function a : X — > Z p , 

2. compute /(A[y]|o-) in tf(n) time, 

3. check whether or not this value is zero and return the result. 

For the time bound note that i/(n) is obviously > \X\. For the error bound if the degree of 
f(A[y]) is d, Corollary | shows the probability of a false zero is < d/p < df(n)/p. □ 

Theorem 9. // the result of f has degree bounded by some polynomial df{n) and f{A[y\\ a ) is 
computable by a straight-line program in tf(n) > df(n) time, then there exists an algorithm that in 
0(tf(n)) time checks for each x £ X whether or not -^ r f(A[y]) is symbolically non-zero over Z p . 
All \X\ returned results are correct with total error probability < df(n)\X\/p. 

Proof. Assume that the Baur-Strassen theorem gives a routine to compute d x f := ^.f(A[y]) for 
all x G X in 0(tf(n)) time. The algorithm works as follows: 

1. uniformly at random choose a substitution function a : X — > Z p , 

2. apply the Baur-Strassen routine to compute d x f\o- for all x G X in 0(tf(n)) time, 

3. for each x £ X check whether or not d x f\ a is zero and return the result. 

Obviously d x f has degree < df(n). So Corollary [3] shows that for each x G X the probability 
of a false zero for d x f is < df(n)/p. Hence the union bound shows the probability of any false zero 
in |X| results is < df(n)\X\/p. □ 

The true goal is to check if a functional value is symbolically non-zero with high probability. 
This is easy to do with some weak assumptions on /. Specifically assume that any value /(a) in 
the range of / is a polynomial whose constant coefficients all have absolute value at most C, for 
some constant C. Also assume tf(n) > n + df(n). 

Now consider the setting of Theorem [8j /(a) is a non-zero polynomial iff it is a non-zero poly- 
nomial over Z p for any prime p > C. In Theorem [8] take p as a prime of size Q(n ■ df(n)). We get 
error probability 0(l/n). The time to find p is easily accounted for by the assumption on tf(n). 
So for / as above, the theorem shows that in 0(tf(n)) time we can check if /(j4[y]) is symbolically 
non-zero with high probability. 

Next consider the setting of Theorem |9| Any value of d x f is a non-zero polynomial iff it is a 
non-zero polynomial over Z p for any prime p > Cd f (n) . In Theorem [9] take p as a prime of size 
G(n| X\ -df(n)). We get error probability 0{l/n). The time is similar. So for / as above the theorem 
shows that in 0(tf(n)) time we can check if J^f(A[y]) is symbolically non-zero with high probability. 
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6 Minimum Weight Perfect Matching 



This section presents an algorithm that, given an undirected graph with integral edge weights in 
[0, W], finds a minimum weight perfect matching in 0{Wn u ) time, assuming such matchings exist. 
The algorithm works in three phases: 

1. The first phase uses algebra to reduce the problem to connected graphs, where each edge be- 
longs to some minimum weight perfect matching (Algorithm [3]) . Moreover for each vertex v, we 
are given the value w(M(v)) - the minimum weight of a matching with exactly 2 unmatched 
vertices, one of which is v (Algorithm [4]) . This phase uses 0{Wn w ) time and succeeds with high 
probability. 

2. The second phase defines a new weight w'(uv) := w(uv) + w(M(u)) + w(M(v)) for each edge 
uv. It performs a simple graph search algorithm on these new edges to obtain a laminar family 
of blossoms, which is the support of some optimum dual solution (Algorithm [5]) . Each blossom 
induces a factor critical graph. This phase is deterministic and uses 0(n 2 ) time. 

3. The last phase uses a maximum cardinality matching algorithm (for unweighted graphs), guided 
by the structure of the blossoms, to obtain a minimum weight perfect matching (Lemma 20). 
This phase uses 0{n u ) time and succeeds with high probability. 

To elaborate on the second phase (which in our opinion is the most interesting), let A be the 



set of distinct values of the weight function w 1 . Section 6.2 proves \A\ = 0{n). For a £ Awe define 
a 'threshold graph' G a = (V,E'), which is an unweighted undirected graph with E = {uv G E : 
w'{uv) < a}. The nontrivial connected components of all the graphs G a constitute the blossoms 
of an optimum dual solution! (A connected component is nontrivial if it has more than 1 and less 
than n vertices.) Our proof of this result hinges on showing there exists a special dual solution 
(called balanced critical dual solution) in which it is easy to find the blossoms (Lemma 28). 

Figure [T] depicts a sample graph and illustrates the steps for obtaining the laminar family of 
blossoms. 




Oi d2 Oi £82 ao Oi a 2 b a bi b 2 c o f-l c-2 



Figure 1: The matching algorithm: The far left shows an example graph. Edges without a label 
weigh 0. The minimum weight of a perfect matching is 3. w(M(v)) is 2 for v G {62, c%, x, y, z} 
and 1 otherwise. The middle figure shows the allowed edges (edge c\z was removed) and new edge 
weights w'{uv) = w(uv) + w(M(u)) + w(M{v)); edges with w'(e) = 2 are drawn zigzag, w'{e) = 3 
are dashed, and w'{e) = 4 are straight. The far right shows the laminar family induced by the 
blossoms found using the threshold graphs of Algorithm [BI 
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In the rest of this section we show how to obtain the set of allowed edges, i.e., edges belonging 



to at least one minimum weight perfect matching and values w(M(u)) (Section 6.1). Next, in Sec- 



tion 6.2 we present the standard LP formulation of the problem. We recall and extend properties 



of a dual solution, in order to prove the correctness of our simple method of obtaining blossoms of 



an optimum dual solution. Finally, in Section 6.3, we gather all the theorems and formally prove 



correctness and bound the running time of our algorithm. 

We would like to note, that the simplicity of our algorithm for contracting the set of blossoms 
from the values w(M(u)) is due to the fact that the hardness is hidden in the proof of the purely 



combinatorial existential lemmas from Section 6.2 and in the algorithm for finding unweighted 
maximum matching problem. 

6.1 Algebraic Tools 

Let us define a symbolic adjacency matrix of the weighted undirected graph G = (V, E, w, W) to 
be the n x n matrix A{G) such that 

{%i,jy w ^ if ij £ E and i < j, 
-Xj ti y w{ij) if ij G E and i > j, 
otherwise, 

where Xij are unique variables corresponding to the edges ij G E of G. Karp, Upfal and Wigder- 
son |19j proved that the smallest degree of y in det(A(G)) is twice the weight of a minimum weight 
perfect matching in G. By using this line of reasoning together with results of Storjohann and 
Baur-Strassen, we show how to obtain the set of edges which appear in at least one minimum 
weight perfect matching. 

Lemma 10. An edge ij G E belongs to some minimum weight perfect matching iff 





det(A(G)) 



J h3 

Proof. By the definition of a determinant we have: 



^0. 



det(I(G)) = Y, S ^(P) II A (G)k, Pk , (2) 
per n k=i 

where T n is the set of n-element permutations. A permutation p defines a multiset of edges 
C p = {{i,pi} : 1 < i < n and i ^ pi}. This edge multiset corresponds to a cycle cover given 
by the cycles of p. Reversing an odd cycle in a permutation does not change its sign, but it changes 
the sign of the monomial corresponding to this permutation in the determinant. Consequently, 
the polynomial det(A(G)) contains only monomials corresponding to even-cycle covers of G. (An 
even-cycle cover has no odd cycles.) Since each even-cycle cover is easily decomposable into two 
perfect matchings, and doubling a matching gives an even-cycle cover, we infer that an edge ij G E 
belongs to some minimum weight perfect matching in G iff it appears in some minimum degree 
monomial in det(A(G)). The lemma follows. □ 

Note that an even-cycle cover corresponding to a matching - i.e., every cycle has length 2 - 
corresponds to a unique permutation in Q. Thus the lemma is true in any field Z p . 
The combination of Lemma 10 with Theorem [9] proves the following corollary. 
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Corollary 11. For a weighted undirected graph G = (V, E, w,W) one can compute the set of edges 
which belong to at least one minimum weight perfect matching in 0{Wn w ) running time, with high 
probability. 



Algorithm 3 Computes the set of allowed edges in the graph G. 



Generate a random substitution a : X — > Z p for a prime p of order 0(n 2 ). 
Compute d = deg* det(A(G)\ a ) using Storjohann's theorem. 

Let d x f be the routine given by the Baur-Strassen theorem to compute the matrix of partial 



derivatives 



d 



term? 



det(i(G)) 



Generate a random substitution a : X 
Compute the matrix 6 = d x f\ a . 
Mark each edge ij, where i < j, as allowed if 5ij 



Z p for a prime p of order 0(n 4 ) 
^0. 



Definition 12 (M(uv), M(u)). For a pair of vertices u, v G V let M(uv) be a minimum weight 
perfect matching in G\ {u, v}, i.e., G with vertices u and v removed. Similarly, for a vertex u let 
M{u) be a minimum weight almost-perfect matching in the graph G \ {u} . 

Note that M(u) always exists, since we assume the given graph G has a perfect matching. In 
contrast M(uv) needn't exist. In that case M(uv) is oo. 

Lemma 13. Let G = (V, E, w,W) be a weighted undirected graph. Then deg* (adj(j4(G))jj) = 
w(M) +w(M(ij)). 

Proof. We have adj(i(G))ij = (-1)*+^ det(i(G) J ''). Equivalently, if we take Z to be the matrix 
obtained from A(G) by zeroing entries of the j'th row and the i'th column and setting the entry 
(j, i) to 1, then &d}(A(G))ij = det(Z), and so 

n 

adj(i(G)) M = ^sgn^J] (3) 

per„ k=i 

The permutation p can be viewed as a set of directed cycles C covering G, where there is a cycle c 
that contains the edge (j, i). (By definition is an edge of C, even when (j, i) ^ E{G). is 
an example. The other edges of C are in E(G).) 

We claim that the terms in this adjoint correspond to even-length-cycle covers that contain a 
cycle c through (J,i). In other words cycle covers that contain an odd cycle make no net contribu- 
tion to In proof let d be an odd cycle in C. First suppose d ^ c. Reversing its direction changes 
the sign of its contribution to ([3]), from antisymmetry of A{G). So such covers do not contribute to 
the adjoint. Second suppose d = c. G has an even number of vertices. Thus C has an even number 
of edges, and c cannot be the unique odd cycle. 

Now take any C contributing to Q . C decomposes into two matchings by taking alternate edges 
from each (even) cycle. One matching, say N, is a perfect matching of G; the other is a perfect 
matching of G \ {i, j}, say N(ij), plus edge (j, i). We conclude that the adjoint of ^ is if N or 
N(ij) does not exist. This proves the case of the lemma when M or M(ij) does not exist. 

Now assume both M and M(ij) exist. The degree of any term in ^ equals the sum of the 
edge weights in C (note that the weight of the edge (j,i) is considered as 0, since Zjj = 1). Take 
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any term that has the smallest degree in ([3]). Its degree is 

degJ(adj(A(G))y) = w(N) + w(N(ij)) > w(M) + w(M(ij)). 

We will complete the proof by showing 

de g ;(adj(i(G)) ij ) < w(M) + w(M(ij)). 

(Obviously the two displayed inequalities show equality holds, so they imply the lemma.) The 
multiset M U M(ij) U (J, i) gives an even- length-cycle cover Cm of G which contains the edge (J, i). 
Wlog the cycle through (j, i) is the only cycle in Cm (since we can make M and M(ij) identical 
on any other even cycle). We claim the monomial corresponding to Cm occurs exactly once in (pjl). 
This claim shows the monomial does not get cancelled, thus proving the desired inequality. 

To prove the claim, a variable x g h appearing in any term comes from the cycle cover edge (g, h) 
or (h,g). This implies that two terms with the same variables are cycle covers that differ only in 
the orientation of some cycles. A cycle of Cm either has length two or contains The former 

does not change when it is reversed, and the latter cannot be reversed. So no other cycle cover 
corresponds to the monomial of Cm- 

Note that since the monomial for Cm occurs exactly once, the lemma holds for any field Z p . □ 

If M or M(ij) does not exist, the expression of the lemma equals oo. The lemma and the following 
extension also hold in any field Z p . 



Corollary 14. The vector of values w(M) + w(M(i)), i 6 V equals deg* (adj(A(G))&) for b a 
vector of n indeterminates b = (pi, 62, . . . , b n ). 



Proof. Lemma 
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; shows w(M)+w(M(i)) = min, deg* [sA}(A(G))ijj . (Note that entries &dj(A(G))ij 

corresponding to the polynomial cause no problem.) The ith component of the vector adj(A(G))b 
is Y2j a dj(A(G))ij6j, and because of the indeterminates bj no terms cancel when the sum is formed. 

Thus deg* (adj(A(G))b\ = w(M) + w(M(i)). □ 

This leads to the following algorithm to compute w(M(u)) for all u £ V(G). Let B = 
{61, . . .,&„}. 



Algorithm 4 Computes the values w(M(i)) for all vertices i in the graph G. 

1: Generate a random substitution o : X U B — > Z p for a prime p of order G(n 3 ). 

2: Compute w(M) = deg*(det(A(G)| (T ))/2 and v = det(i(G)| (T )(i- 1 (G)6)| (7 using Storjohann's 

theorem, where b = (61, ... , b n ). 
3: For each % £ V set w(M(i)) = deg*(wj) - w(M). 
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To see this algorithm is correct Statement 2 computes v = ( adj(A(G))6) | CT . So Corollary 
shows the algorithm is correct if there are no false zeroes. A rational expression (like those in 
(A _1 (G)6)| CT ) is zero if and only if its numerator is zero and its denominator is nonzero. So we can 
apply the Schwartz-Zippel Lemma to show the final products have no false zeroes (in their lowest 
order term). So each deg*(t>j) is computed correctly with high probability. 
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Regarding efficiency consider the n multiplications of degree nW polynomials done to form v in 
Statement 2. We only use the lowest degree term of each product (Statement 3). That term comes 
from the lowest degree term in det(A(G)| CT and the lowest degree term in the numerator and the 
denominator of (A~ 1 (G)b)\ a . So we can find the smallest degree of y that corresponds to deg*(vi) 
using 0(1) additions and subtractions, without multiplying polynomials. 

We conclude: 

Corollary 15. Algorithm^ computes the values w(M(u)), for all u £ V(G), in OiWn UJ ) time, 
with high probability. 



6.2 Properties of the Dual 

We move on to the linear programming formulation of the minimum weight perfect matching prob- 
lem given by Edmonds [5]. An odd set has odd cardinality; Q denotes the collection of odd subsets 
of V of cardinality > 3. 



min w(e)x e 

x(8(v)) = 1, for all v £ V 

x(5(U)) > 1, for all U £ (4) 
x e > 0, for e £ E 

The variables x e indicate when an edge is included in the solution. Here, 5(U) denotes all edges 
uv G E having \{u, v} f]U\ = l. We write 5(u) for 5({u}) and x(F) for ^2 ee p x e- 
The dual problem has variables ir v for each vertex v and ttu for each odd set U: 



max > tt v 



tt u + tt v + ttu < w{uv) for all uv G E (5) 

uen, uves(u) 

ttu > for all U e 

We say that an edge e = uv is tight with respect to a dual it if equality holds in . A laminar 
family is a set system where each pair of sets is either disjoint or one set contains the other. More- 
over, a graph is factor critical if after removing each vertex the graph has a perfect matching. We 
use existence of the following dual: 

Lemma 16 (Edmonds '65 [5]). There exists an optimal dual solution tt : V U Q — > M, such that: 

1. the set system {U £ 0, : it\j > 0} forms a laminar family, 

2. for each U £ O with nu > 0, the graph G[U] with each set of {S £ £1 : S C U, tts > 0} contracted 
is factor critical. 

Definition 17 (critical dual, blossom). An optimum dual solution satisfying the conditions from 
Lemma 16_ is a critical dual solution. A set U £ Q such that ttjj > is a blossom w.r.t. tt. 

Blossoms of critical dual solutions have the following useful property (note the lemma below is 
weight-oblivious and the only input given to the algorithm is an undirected unweighted graph, the 
family of blossoms, and v). 
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Lemma 18. Consider any critical dual solution and let U £ H be an arbitrary blossom. For any 
vertex v £ U there exists a perfect matching M(U,v) in G[U \ {v}], such that for each blossom 
Uq C U, \M(U,v) n 8(Uq)\ is if v G Uq and 1 if v G" Uq. Furthermore, given the family of all 
blossoms and v, one can find such a matching in 0{\U\ U ) running time, with high probability. 

Proof. Let B be the set of blossoms of ir properly contained in U {B might be empty); moreover 
let £> max be the set of inclusionwise maximal sets in B. Let G' be the graph G[U]/B max and let v' 
be a vertex of G' corresponding to v. (Here we use the contraction operator - if S is a family of 
disjoint vertex sets, G/S denotes the graph G with each set of S contracted to a single vertex.) 

Initially let M(U, v) C E be a perfect matching in G'\v'. It exists since G' is factor critical. For 
each blossom Uq G B m ax, recursively find a perfect matching Mq in the graph G[Uq \ x], where x is 
the single vertex of the intersection of Uq and V(M(U, v)) U {v}. Add the edges of Mq to M(U, v). 

By construction the final matching M(U, v) satisfies conditions from the lemma. For the time 
bound note that the laminarity of B implies the total number of vertices in all graphs constructed 
by the above procedure is 0(n). The algorithms of [2'6\ E?o] IT4"] find a perfect matching on an 
arbitrary graph of n vertices in time 0{n w ), with high probability. Hence our recursive procedure 
runs in total time 0(| f/| w ). □ 

Complementary slackness gives the following observation. 

Observation 19. For any optimum dual solution: 

(a) a set U E Q with ttjj > has exactly one edge of 5(U) in any minimum weight perfect matching; 

(b) an edge belonging to any minimum weight perfect matching is tight. 

Lemma 20. Given a weighted undirected graph G = (V,E,w,W) where each edge is allowed, 
and the set of blossoms B of some critical dual solution, one can find a minimum weight perfect 
matching in 0{n w ) time, with high probability. 

Proof. Let £> max Q £> be the set of inclusionwise maximal blossoms. Let graph G' = G/B max . G' 
has a perfect matching Mq (since Observation [l9|a) shows any minimum weight perfect matching 
in G contains a subset of edges forming a perfect matching in G'). 

We extend Mq to a perfect matching in G by considering the blossoms of U G £> max one by 
one. Let vjj be the unique vertex of U that is matched by Mq. Add the matching M(JJ,vjj) of 



Lemma [18] to Mq. The final set Mq is a perfect matching for G. 

since the total 
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The time for this procedure is 0{n u ). This follows exactly as in Lemma 
number of vertices in all graphs considered is 0{n). 

Finally, we prove that M is a minimum weight perfect matching in G by showing that for 
each blossom U G B the set M contains exactly one edge of 6(U). Consider a maximal blossom 
U G i3 max . When finding a perfect matching in G' we have added exactly a single edge of S(U) to 
the set M. Moreover each edge of M(U',Vu'), for any maximal blossom U' G £> max , is contained 
in U', and therefore the set M contains exactly one edge of 8(U). Now let us consider a blos- 
som U G B \ B max , and let U' G B ma be a maximal blossom containing U . If in the first phase 



(finding a perfect matching in G 1 ) we added no edge from 5(U) to the set M, then by Lemma 18 
in the set M(U' %/) there is exactly one edge of 5(U), whereas for each other maximal blossom 
U" G £> max , U" / U', in the set M(U" ,vjj») there is no edge of S(U). If, however, in the first phase 
we added an edge from 5(U) to the set M, then by the choice of v'jj, which is the endpoint of the 
edge of S(U) n M, no other edge of 5(U) is added to the set M. 
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Since all the edges are allowed, by Observation 19 b) the sum of values of edges of M is equal 
to the cost of the critical dual solution, which proves that M is a minimum weight perfect matching 
in the graph G. Note that our algorithm does not need exact values of the dual nor even weights 
of edges, since its input is only graph G and set B. □ 

A critical dual solution gives rise to a weighted tree in a natural way: 

Definition 21 (dual tree). Let ir : £lL)V — > M. be a critical dual solution, with B the set of its 
blossoms. The dual tree T(ir) is a rooted tree on nodes {V} U BU V , where V is the root, vertices 
ofV are leaves, blossoms of B are internal nodes and the parent-child relation in T(ir) is naturally 
defined inclusionwise. The weight of the edge from a node t G BL) V to its parent is irt- The height 
of the tree H(T(ir)) is the weight of a longest path from the root to some leaf. 

In this definition note that the last edge of a path defining H(T(jr)) may have negative length. 
For a tree T with weighted edges and two nodes u, v, distx(u,v) denotes the weight of the path 
between u and v. The following simple lemma provides a basic tool. 

Lemma 22. If ir is a critical dual solution for a weighted graph G = (V, E,w,W), any allowed 
edge uv satisfies w(uv) = dist T ^(u,v). 



Proof. Since uv is tight (Observation 19 'b)), w(uv) = ir u + ir v + Ylueti uve8(u) ^u- The right-hand 
side gives dist^^) { u -> v ) f° r t wo reasons: The edges of T(n) incident to leaves are weighted with the 
singleton values of ir. A blossom B of it contains exactly one endpoint of the edge uv if and only 
if the path between u and v in T(jr) contains the edge between B and its parent. □ 



The next steps of our development (Lemmas 25 -27) can be derived using an appropriate ver- 



sion of Edmonds' weighted matching algorithm (e.g., |29J). Here we will use a structural approach, 
based on the following properties of allowed edges given by Lovasz and Plummer. 

Lemma 23 ([21). Lemma 5.2.1 and Theorem 5.2.2). Let G = (V,E) be an undirected connected 
graph where each edge belongs to some perfect matching. Define a binary relation R C V x V by 
(u, v) £ R if and only if G\ {u, v} has no perfect matching. Then 

• R is an equivalence relation; 

• each equivalence class of R is an independent set; 

• for each equivalence class S of R, the graph G\S has exactly \S\ connected components, each of 
which is factor critical. 

We will use a special type of critical dual solution that we call " balanced" . 

Definition 24 (balanced critical dual). Let 7r : fiU 7 - > M be a critical dual solution, and 
let G' be the graph G with each blossom of tt contracted, ir is a balanced critical dual solution if 
there are two distinct vertices u,v E V such that distT( n )(u,V) = dist T ^(v,V) = H(T(ir)) and 
further, G' \ {u',v'} has a perfect matching for u',v' the (distinct) vertices of G' corresponding to 
u,v, respectively. 

Before proving that balanced critical dual solutions exist, we give a lemma showing why they 
are useful. In particular they show how the M(y) values relate to T(ir). Let M(G) be a minimum 
weight perfect matching in G. 
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Lemma 25. Let G = (V, E, w, W) be an undirected connected graph with every edge in some mini- 
mum weight perfect matching. Let ir be a balanced critical dual solution for G. For any vertex z G V , 
a minimum weight almost perfect matching in G\z weighs w(M(G)) — H(T(jr)) — distxu)[z, V). 

Proof. Any almost perfect matching in G \ z weighs at least w(M(G)) — H(T(ir)) — dist T ^(z, V). 
In proof let M\ be an arbitrary perfect matching in G \ {x, z} for any x £ V. For any blossom U 
of vr such that x,z U, \Mi n 5{U)\ > 1. Together with ^ this gives w(M{) > Y^weV-x^^ + 
YxMV n U- The right-hand side equals w(M(G))— ir x — ir z — J2{ x ,z}r\U^d ^Ui strong duality. Since 
every iru is nonnegative this quantity is at least w(M(G)) — distxu) (z, V) — distj>u) (%, V). The defi- 
nition of H(T(it)) shows the last quantity is at least w{M{G)) — dist T r n \(z, V) — H(T(ir)) as desired. 

We complete the proof by constructing an almost perfect matching in G\z of weight w{M{G)) — 
H(T{n)) — distx^)(z,V). Take G' ,u,v,u' ,v' as in Definition 24 Moreover let z' be the vertex of 
G' corresponding to z. G' is connected, with every edge in a perfect matching, so it satisfies the 



hypothesis of Lemma 23 Definition 24 shows that u' Ifiv 1 . So z' is not equivalent to at least of u 
and v . W.l.o.g. assume that u' Ljlz 1 . Thus G' \ {u' , z'} has a perfect matching Mq. 

Next, consider each inclusionwise maximal blossom U of tt one by one. Let x E U be the unique 
vertex of U in the set V(Mq) U {u,z}. Add to Mq the edges of the matching M(U,x) guaranteed 



by Lemma 18 



Clearly Mq is a perfect matching in G \ {u,z}. For each blossom U of n, \Mq n S(U)\ is 1 
if u, z U, and if u or z belongs to U. Blossoms of the latter type are those in the path 
from u to V or z to V in T(n). These two paths have disjoint edge sets, since v! ^ z' . We get 
w{Mq) = w(M(G)) — distT( 7T )(u,V) — dist T ^(z,V), since every edge of Mq is allowed, i.e., tight. 
Since dist T ^(u, V) = H(T(tt)) this is the desired weight. □ 

We prove that balanced critical duals exist in two steps. The first step shows a simpler property 
for critical duals actually makes them balanced. The second step shows duals with this property 
exist. 

Lemma 26. Let G = (V, E, w, W) be an undirected connected graph with every edge in some 
minimum weight perfect matching. A critical dual ttq is balanced if it has minimum height (i.e., 
H(T(ttq)) is no larger than the height of any other critical dual). 

Proof. Assume for the purpose of contradiction that ttq is not a balanced critical dual. For any 
vertex v S V let h v denote its height in ttq, h v = dist T ^ ){v, V). Let u be the vertex of G with the 
greatest height h u . Let G' be the graph G with inclusionwise maximal blossoms of ttq contracted. 



Let R be the equivalence relation of Lemma 23 for G', and S%, . . . , Sk its equivalence classes. Let 
u belong to vertex u' of G' and let u' £ Si. 

We will define a dual function ~k\. An element of S\ is either a maximal blossom of ttq or a 
vertex of V not in any blossom; let Sj, 1 < i < be the zth of these blossoms and vertices. 



Lemma 23 shows G' \ Si has |«Si| connected components; let Bi, 1 < i < | S'i | , be the set of vertices 
of G contracted onto the i-th connected component of G'\S\. Define tv\ : f2 U V — > IR to be identical 
to ttq except 

J vro(x) — e x = Si, 1 < i < \S\\ 
1 ttq(x) + e x = Bi, 1 < i < \S±\. 

(Note that if Bi consists of more than one vertex in G' then we are creating a new blossom.) Let 
e be any positive real no larger than the smallest value of 7ro(si) for a blossom Sj. This ensures tt\ 
is nonnegative on blossoms. 



TTi(x) 
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Figure 2: Graph C and the modification of the duals. 



Let us verify that %\ is a critical dual. First, observe that each edge of G remains tight in tt\: 
Nothing changes for an edge that has both ends in the same set of S± or some B{. The remaining 
possibility is an edge between Si and some Bi (no edge joins 2 Bi sets or 2 sets of S%, the latter 
by independence of Si). For such edges we have added and subtracted e in the left-hand side of[5j 



so it remains tight. Next observe that the blossoms of tti form a laminar family. Lemma 23 shows 
the sets Bi induce factor critical graphs. Finally m is an optimum dual, since its objective as ttq. 
Thus it i is a critical dual. 

Taking e small enough makes tti a critical dual with smaller height than ttq, the desired contra- 
diction. To see this take any vertex v E V, and let v' be the vertex of G' that v is contracted onto. 



If v' G Si, the height of v decreases as long as e is positive. Suppose v' Si- Lemma 23 shows 
7To would be balanced if h v = h u . Thus h v < h u . Choose e small enough so that every such v has 
dist T ^^{v) = dist T i m \{v) + e < h u — e. Thus m has smaller height than ttq. □ 

Lemma 27. Let G = (V,E,w,W) be an undirected connected graph with every edge in some min- 
imum weight perfect matching. There is a critical dual ttq that has the smallest height H(T(ttq)). 

Proof. Lemma [16] shows a critical dual n exists. There are a finite number of laminar families on 
V, i.e., a finite number of trees T(ir). So it suffices to show that there is a smallest height among 
all critical duals tt with the same tree T = T(tt). 

We begin by showing that for every blossom U, there is a unique value for tt x , where x is any 
vertex of U or any blossom properly contained in U. We argue inductively, so assume this holds 
for every blossom properly contained in U. For any edge uv we break the left-hand side of ^ into 
the contributions from u and from v, by defining 

UeQ, u={u,v}nu 



l v,u- 



and symmetrically for % VjU . So the left-hand side of ([5]) is tt U:V + ir v 

Take any edge uv joining two vertices u,v £ U. uv is on an odd cycle C contained in U. (U is fac- 
tor critical, so let M u (M v ) be a perfect matching on U—u (U—v) respectively. The symmetric differ- 
ence M U ®M V contains an even-length path from u to v.) Each edge of C is tight. So for every edge 
xy in C, the values of ir x ,y and ir yiX are uniquely determined. If ir x , y does not have any contributions 
from blossoms properly contained in U then tt x = tt x ^ has been uniquely determined. If n Xt y has 
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a contribution ttw from a blossom W that is a maximal blossom properly contained in U then ttw 
has been uniquely determined. This follows since the other tt values contributing to Tr x ,y have been 
determined by induction. (Note that ttw has also been uniquely determined from the other edge of 
C n 5(W).) If neither of these conditions apply to ir Xl y then all its -/r- values have been determined 
by induction. Since any vertex u E U is on an edge uv in U, this completes the inductive argument. 

Next consider any edge uv not contained in a blossom of T. The previous argument shows 
exactly one term in the quantity 7r Ujl , is still undetermined. If uv is in an odd cycle C the previous 
argument shows that term is uniquely determined. Contract all such odd cycles as well as all blos- 
soms of T. We get a bipartite graph G' . It contains at least one edge. Let S be a spanning tree of 
G'. Choose a value po for the unknown term p at the root of S, that comes from a valid critical dual 
for T. Suppose we increase p. If this increases H{T), every value of p larger than pq gives larger 
height. Suppose this decreases H(T). All the other unknown tt- values are uniquely determined 
from tightness of the edges of S. Also every edge of G' not in S remains tight by bipartiteness. 
There is a maximum value p such that every value p > p either makes the 7r-values invalid (because 
some ttu, U G becomes negative) or increases the height (since p contributes to the height of the 
root vertex). Similarly there is a minimum value p for p. We conclude there is a unique smallest 
height for a critical dual for T - it occurs when p is equal to either p or p. □ 

As already mentioned, the last two lemmas show any undirected connected graph G = (V, E, w, W) 
with all edges allowed has a balanced critical dual. We can now reach our final goal. 

Lemma 28. Let G = (V, E, w, W) be a weighted undirected connected graph where every edge is 
allowed. Given all values w(M(v)) for v E V, the blossoms of a balanced critical dual solution can 
be found in 0(n 2 ) time. 



Proof. Let tt be a balanced critical dual solution. By Lemma 25 for each leaf node v € V of T{tt), 
w(M(v)) = w(M(G)) - H(T(it)) - dist TM (z,V). Define new edge weights w' : E -> Z as w'{uv) = 
w(M{u)) +w{M(y)) +w(uv). Consider any uv £ E. Since uv is tight, w(uv) = dist T ^(u,v). 
Define a quantity c that is independent of uv, c = 2{w{M{G)) — H(T(tt))). Then 

w'(uv) = w(uv) + w(M(u)) + w(M(v)) 

= dist T ( w) (u,v) + 2w(M(G)) - 2H(T(tt)) - dist T{7r) (u,V) - dist TM (v,V) 

= 2{w{M(G)) - H{T(tt))) - 2dist T(w) (lca(u,v),V) 

= c — 2distxu)(lca(u,v),V) . (6) 

Let B = lca(u,v). So B is the inclusionwise minimal blossom of tt containing both u and v, or 
if no such blossom exists, B is the root V of the tree T(tt). For any edge uv let B uv C V be the 
the set of vertices reachable from u or v by a path of edges e satisfying w'{e) < w'(uv). 

Claim. For any edge uv, B uv = B. 

Proof of Claim. Let F C E be the set of edges of a spanning tree of G[B] (G[B] is connected since 
either B = V or G[B] is factor critical). Since any edge ab £ F is contained in B, the node lca(a, b) 
descends from lca(u,v) in T(tt). Thus the path from lca(a,b) to lca(u,v) in T has nonnegative 
weight. This implies w'(ab) < w'(uv) by Q. Thus B C B uv . 

For the opposite inclusion, consider any edge ab with a £ B and w'{ab) < w'(uv). Since every 
blossom has a strictly positive 7r-value, Q implies b G B. Now an easy induction shows any path 
from u or v, with every edge e having w'{e) < w'(uv), has every vertex in B. Thus B uv C B. 
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Any blossom B of ir has an edge uv with B the minimal blossom containing u and v (by 1am- 
inarity and connectedness of B). So the claim of the lemma amounts to constructing all the sets 
B uv . This is done in 0(n 2 ) time by Algorithm pi below. □ 



Algorithm 5 Given all the values w(M(u)), finds the blossoms of a balanced critical dual in the 
graph G where all edges are allowed. 

1: For each edge uv set w'{uv) = w(uv) + w(M(u)) + w(M(v)). 

2: Let A be the set of all different values w'(uv). Let B = 0. 

3: for each a G A, in increasing order, do 

4: Let C be the set of connected components of the graph (V, {uv : uv £ E, w'(uv) < a}). 
5: Add the nontrivial components of C to B. 
6: end for 
7: return B. 



6.3 The Final Algorithm 

Theorem 29. Let G = (V, E,w, W) be a weighted undirected graph containing a perfect matching. 
A minimum weight perfect matching in G can be computed 0{Wn u ) time, with high probability. 



Proof. First, using Corollary we can remove all the edges of G which are not allowed. Clearly, 
we can consider each connected component of G separately, hence w.l.o.g. we assume that G is 



connected. Next, compute all the values w(M(u)) for each u £ V using Corollary 15 Having all 



the values w(M(u)) by Lemma 28 we can find the set of blossoms B of a balanced critical dual 



solution and consequently by Lemma 20 we can find a minimum weight perfect matching in G. □ 



The full version of this paper shows how the matching algorithm can be made Las Vegas. 

In some applications the second smallest perfect matching is of interest (e.g., Section [9]). Its 
weight is easily found, as follows. As discussed in the proof of Lemma [T3] , the terms in the de- 
terminant of det(A(G)) correspond to even-cycle covers in the graph G. Each such cycle can be 
decomposed into two perfect matchings. As we already observed the smallest degree term in y 
in det(A(G)) corresponds to taking twice the minimum weight perfect matching in G. The next 
smallest term gives the following. 

Corollary 30. Let G = (V, E, w, W) be a weighted undirected graph. The degree in y of a second 
smallest monomial o/det(A(G)) is equal to the weight of a minimum weight perfect matching M* 
plus the weight of a second smallest perfect matching M' . In particular, the weight of a second 
smallest perfect matching can be found in 0{Wn u ) time, with high probability. 



Proof. Letting r(M*) denote the term corresponding to M* in det(A(G)), 

w(M') = term* fdet(A(G)) - r(M*) 



□ 



7 Diameter and Radius 

In this section we consider the problem of computing the diameter and radius of a directed graph 
without negative weight cycles. By bidirecting the edges this result can be applied to undirected 
graphs with nonnegative edge weights. 
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We start with definitions of the quantities of interest, plus equivalent definitions that we use to 
compute these quantities. To motivate the latter note that it seems difficult to compute a given set 
of distances S directly. Instead we show how to check if, for an arbitrary value c, all the distances 
in S are < c. For a graph G and a vertex i, 

eccentricity (i) = max{disic(i, j) : j £ V} = min{c : (Vj)(c > dist(i,j))}, 
radius(G) = mm{eccentricty(i) : i £ V} = min{c : (3i)(Vj)(c > dist(i, j))}, 

diameter(G) = m&x{eccentricty(i) : i E V} = min{c : (Vi, j)(c > dist(i,j))}. 



We use the following theorem proven in [27jj 
Lemma 31. Let G be a directed weighted graph without negative weight cycles. The weight of the 
shortest path in G from i to j is given by dista(i,j) = deg* ^adj ^A0) + 1 ^ V Moreover, all 

non-zero terms in det (a(G*)) are non-zero over any finite field Z p . 

In order to be able to use the above lemma we first need to apply the following observation. 
Corollary 32. Let c be arbitrary number from [— nW, . . . , nW] . There exists d < c such that 

/ 0. This 



term^ 

L x 7 J J 

continues to hold in any finite field 



adj (A0) . + if and only if term^ adj (A0) +*)..- (EST V' 



Proof. Multiplication of a polynomial by Ei=cT y % means that we add all lower degree terms to a 
term of a given degree. Hence, if some degree term d was non-zero then all higher degree terms 
become non-zero (assuming no terms get cancelled). Since < c — d < 2nW a degree d term creates 
a corresponding degree c term in the product. 

Finally observe that no term drops out because of cancellations^} In the matrix A(G) + L, every 
nonzero entry has the form xy w ^ , where x is the indeterminate 7^ y (or 1 in diagonal entries) and 
w is a function. So every term in + 1) has the exponent of y functionally dependent on 

the remaining variables. This holds for an entry of the adjoint too. Thus a term (Yl x )y c created 
in the multiplication of the corollary comes from exactly one term in the adjoint (i.e., 

d = ^2w(x)). So there are no cancellations, in ordinary arithmetic or in Z p . □ 

To find the diameter, we use this Corollary to perform a binary search for the lowest c such 
that for all i,j G V, 

/2nW 

adj (!(£) + /] 



term^ 




Clearly this c is the diameter. Similarly a binary search for the lowest c where the displayed 
condition holds for some i with every j gives the radius. 

adj (i(^) + l). . • fe^V 



The main problem left is how to check whether term^ adj ( A(G) + / ) • ( E«=o U % I ^ 0- 

We will show how to obtain adj as a partial derivative of det. We define Z to be a fully symbolic 
matrix of size nxnas Zjj = Zij, where Zij are unique variables for all i, j G [1, . . . , n]. We define 
a z to be an evaluation that assigns to all Zij. Now we are ready to prove the following lemma. 



3 In |27j the graph was defined to contain self-loops whereas here we add self-loops in the equation by taking 
A0) + I. 

4 We shall apply this principle to other matrices. 
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Lemma 33. 



term 



f2nW 



adj(i^)+i). .(xy 



. i=0 



det (A0) + I + Z) 



8=0 



Proof. Observe that for the fully symbolic n x n matrix Z and any n x n matrix ^4 not involving 
any variable Zij, 







dz 



det(i + Z)| CTz = adj(i) 



(7) 



J: I 



Thus 

9 d 
a term 



term! 



u 



I znvv 

det (A0) + 



adj (A0) + I + Z 



/2nW \ ' 






>(£'). 




= term^ 









d 



dz 



det (A0) +I + Z 



3,1 



/2nW N 
A i=Q / 



/2nVK 
\ i=0 



term, 



/2nW 



adj(l0) + l).. ^ 



\ i=0 



□ 



Joining the above results together with Theorem [9] and Corollary 32 we get the following algo- 
rithm. 

Algorithm 6 Checks whether diameter of the directed graph G is < c. 
1: Let d z f be the routine given by the Baur-Strassen theorem to compute the matrix of partial 
derivatives gj. term c + nW fegf y*) det ((10) + J + Z)y w \ ] . 

2: Generate a random substitution ct : X — >• Z p for a prime p of order 0(n 4 ). Extend it to 

cr : X U Z — ^ Z p by setting <j|z = a z . 
3: Compute the matrix 5 = d z f\ a . 
4: Return true if 5ij is non-zero for all i,j G V. 



Note that the polynomial multiplication in Step 1 need only compute the coefficient of y c+nW 
and so only uses time 0(nW). Using binary search with the above algorithm, and using a similar 
algorithm for radius, we obtain: 

Theorem 34. Let ~G* = (V,E,w,W) be a weighted directed graph without negative weight cycles. 
The diameter and radius of G can be computed in 0(Wn ul ) time, with high probability. 

8 Shortest Cycles in Undirected Graphs 

Let G = (V, E, w, W) a weighted undirected graph. In this section we consider only the problem 
of computing shortest cycles when all edges have non- negative weights. The more general case is 
solved in the next section, and uses the ideas introduced here. Let us define a symbolic adjacency 
matrix of the weighted undirected graph G to be the be the symbolic matrix polynomial A(G) 
equal to A(G) , where G is the bidirection of G. 
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Lemma 35. Let G be an undirected weighted graph with no negative edges, and d E [0,nW]. Some 



uv E E has 



' rl-r 



term 



7^ iff some cycle packing of weight 



det(A{G) + I) — 1 

d contains a component that is an oriented cycle through uv in G. Moreover, for p > 2, all non-zero 
terms in the above expression are non-zero over the finite field Z p . 

Proof. As argued in the proof of Lemma |4]det( A(G) +1) — 1 contains only terms that correspond to 
cycle packings. Moreover, the degree of each non-zero term is equal to the total weight of the cycles 
in the packing. However, because A{G) was constructed using bidirected graph there might be terms 
containing both antiparallel edges, that correspond to bidirected edges in the undirected graph. 

Next we show that a cycle packing C in G contributes to the expression of the lemma iff it con- 
tains exactly one of the variables x U)V and x V;U , and hence contains a simple cycle passing through 
uv. Moreover we show that in such a case the contribution of C is a product of variables corre- 
sponding to C and hence the contribution of C is not cancelled out by a different cycle packing. If 
C does not contain 

%u,v or x v u , clearly it has zero contribution. Xhis leaves two possibilities: 



Case 1. C contains a cycle u,v,u Since 
0, C's term makes no contribution. 



x 



d 



U.V ftr 



— X 



d 



v > u dx v . 



•Eu,v%v : u 



Xu.v^i: 



Case 2. C contains a simple cycle C containing uv The corresponding term contains exactly 



one of x uv and x vu say x uv . We have 



x uv . Hence, the derivative 



for this term is nonzero and is equal to the sign of permutation multiplied by the product of the 
variables of the oriented edges of C. □ 

Similarly as in undirected graphs we say that edge e is allowed if and only if it belongs to some 
simple shortest cycle C in G. The above proof actually gives us a way to find allowed edges as well. 

Corollary 36. Let d be smallest number in [l,nW] such that there exists an edge uv E E such 



that 



d 



only if when 



x 



d 

v,u 

d 



term! 



det(A(G) + /) — 1 0. Then an edge uv E E is allowed if and 



X 



d 



terrn^ 



det(i(G) + /) - 1 



7^0. 



Computing the Weight of the Shortest Cycle Using Lemma 35 we will devise an algorithm 



that will be able to check whether there exists a simple cycle in G of length shorter or equal to c. 
In order to do it we need the observation similar to Corollary |32| 



Corollary 37. Let c be arbitrary number from [l,nW]. There exists d < c such that 

d d 



J, ■,. 



term 



det(A(G) + /)-l 



^0 



if and only if 







d 



terniy 



(det(i(G) +/)-!)•(£ 



o v 



Using the above observation we can construct the following algorithm. 
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Algorithm 7 Checks whether the shortest cycle in undirected graph G has weight < c. 
1: Let d x f be the routine given by the Baur-Strassen theorem to compute the matrix of partial 



derivatives 



term? 



det(i(G) + I) - 1 



Generate a random substitution a : X — > Z p for a prime p of order 0(n 4 ). 
Compute the matrix 5 = d x f\ a . 

Compute the matrix 5' with 5' u v = [x U:V 5 U:V — x v>u 6 VtU \ | . 
Return true if 8' has a non-zero entry. 



The correctness of the above algorithm is implied by both Theorem [8] and Theorem [9} Using 
binary search with it we obtain. 

Theorem 38. Let G = (V,E,w,W) be a weighted undirected graph without negative weight edges. 
The weight of the shortest simple cycle in G can be computed in 0(Wn u ) time, with high probability. 

Finding the Shortest Cycle After showing how to compute the shortest cycle length it re- 
mains to show how to find the cycle itself. We essentially can use the same approach as we used 
for directed graphs in Section |4j 



Algorithm 8 Computes the shortest cycle in undirected graph G. 



Let c* be the weight of the shortest cycle computed using Theorem 38 
Let 5' be the matrix computed by Algorithm [7] for c = c* . 
Take any edge uv such that 5 UV ^ 0. 

Compute the shortest path p VjU from v to u in G \ {uv} using Dijkstra 
Return the cycle formed by uv and p v u - 



Theorem 39. Let G = (V,E,w,W) be a weighted undirected graph without negative weight edges. 
The shortest simple cycle in G can be found in 0{Wn UJ ) time, with high probability. 

9 Undirected Graphs with Negative Weights 

This section gives algorithms for shortest cycle and diameter in undirected graphs with possibly 
negative edges but no negative weight cycles. To accomplish this we need to combine the ideas from 
Sections [7] and [8] with our results for weighted matching. We will recast the results for the diameter 
and shortest cycles into the language of matchings. Hence (unlike Section [8]) throughout this section 



the symbolic adjacency matrix A(G) of an undirected graph G is defined as in Section 6.1 



Diameter Let G = (V,E,w,W) be an undirected graph with negative weights allowed, and let 
E~ be the set of edges with negative weights. We will define a graph G that models paths in G by 
almost perfect matchings. We believe the construction is essentially due to Edmonds [6]. 
Define the split graph G = (V, E, w, W) where 

V = {vx, v 2 : v e V} U {ei, e 2 : e G E~}, 

E = {V1V2 ■ V G V} U {uxV2, U2VX, UxVi,U2V2 ■ uv £ E \ E~} 

U {uiei,it2ei,eie 2 ,fie 2 ,f2e2 : e = uv e E~,u < v}, 
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w(uv) if uv £ E \ E~ , 
w(uiVj) = I w(e) if Ui = ei and Vj ^ e2 and e£i?" 
otherwise. 



&1 ^2 &1 &2 




Figure 3: An undirected graph G and its graph G. In G zigzag edges weigh —1, dashed edges 
weigh 1 and the remaining edges weigh 0. Vertices corresponding to negative edges of G are white 
squares. The far right shows a matching M(a2C\) of weight —2, which corresponds to a shortest 
path between a and c. 



Note how a length-two path in G, say a, b, c with w(ab) > > w(bc) and e = 6c, corresponds to a 
matching in G such as ai&i, 6261, &2Ci, having the same total weight. Figjijgives a complete example. 

An important property is that we can assume n = \ V\ < An. This follows since we can assume 
\E~ \ < n, as otherwise the set of negative edges contains a cycle. 

We will consider minimum weight perfect matchings in G. To use our algebraic tools we should 
eliminate negative weights by setting w'(e) := w(e) + W. Obviously this increases the weight of 
all perfect matchings by hW/2 and so doesn't change the minimum perfect matching. But to keep 
things simple in the following we keep G as defined above, with w possibly negative. 

The following observation is essentially given in pQ in Chapter 12.7 (for a larger version of our 
graph). 

Lemma 40. Let u,v E V, let M be the minimum weight perfect matching, and let M{u2V\) be 
the minimum weight almost perfect matching in G that does not match v\ nor u-i- If G does not 
contain negative weight cycles then w{M) = and the shortest path weight from u to v in G is 
equal to w(M(u2Vi)). 

Note also that it is easy to detect a negative cycle in G - it corresponds to a perfect matching 
in G with negative weight. 



By Lemma 13 we know that deg* (&d}(A(G)) U2tVl ) = w(M) + w(M(u2Vi)). Thus distc{u,v) 



deg* (&dj(A(G)) U2%Vl ), i-e., just as in Lemma 31, adj(vl(G)) encodes the distances in G. Now we 
proceed exactly as in Algorithm [61 
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Algorithm 9 Checks whether diameter of the undirected graph G with negative weights is < c. 



Let d z f be the routine given by the Baur-Strassen theorem to compute the matrix of partial 

^to V) det ((1(G) + Z)y w " 



derivatives 



<->.-. 



2^1 » 



,m,d£ V. 



Generate a random substitution a 
er : X U Z — )• Z p by setting <j\z = o~z- 
Compute the matrix 5 = d z f\ a . 
Return true if S U2jVl is non-zero for all u,v £ V. 



X — > Z p for a prime p of order 0(n ). Extend it to 



Note that ^4(G) + Z is not skew-symmetric but we still get the adjoint by ([7]). Similarly to 
check if the radius is < c, Step 4 returns true if some row of 5 consists entirely of nonzeroes. Again 
using binary search we obtain. 

Theorem 41. Let G = (V, E, w, W) be a weighted undirected graph without negative weight cycles. 
The diameter and radius of G can be computed in OiWn^) time, with high probability. 



Shortest Cycles Recalling Corollary |30[ it might appear that a second smallest perfect matching 
in G corresponds to a shortest cycle in G. But this is not true, because of cycles of length two 
(e.g., uiVi,V2U2 or U2Vi,V2U%). These can be handled as in Section[8j by antisymmetric derivatives 
For the next lemma note that in the absence of negative cycles, any cycle 



contains an edge of nonnegative weight. 

Lemma 42. Let G be an undirected weighted graph with no negative cycle. For any edge uv E 



E\E , a shortest cycle through uv weighs term* 
continues to hold in any field Z p , p > 2. 



„ d _ d 



det(A(G))). This 



Proof. Let 5 =- ^ Ul) v 2m — 
ables x, y and any integers i,j, 



det(_A(G)) ) . Observe that in general for any vari- 



8__ d_ 
dx ^ dy 



x l yi = (i- j)x l y 3 . 



(8) 



Hence the terms in 5 are a subset of those in det(A(G)), i.e., the effect of the differentiation operator 
is just to change the multiplicity of some terms, perhaps zeroing them or causing other cancellations. 

Let C be a shortest cycle through uv. It gives an even-cycle cover M* U N in G with weight 
w(C), where M* is the minimum perfect matching of G, M* = {v±V2, e\e2 '■ v S V, e E E~}, and N 
is the perfect matching of G containing u\V2 but not U2V1, plus representatives of the other edges 
of C, plus edges X1X2 for vertices or negative edges x ^ C. Let r be the monomial corresponding to 
this cover (e.g., r has the term y w ( c >). Q (applied to x^ ^x^ ^) shows r is a term contributing 
to 5. In fact r is the only such term in S involving its monomial. This again follows from ^ and 
the preliminary observation (it is easy to see r is the unique edge cover for its monomial, i.e., none 
of its cycles can be reversed). We conclude term*(5) < w(C). 

We complete the proof by showing term* (5) > w(C). First observe that a perfect matching M 
on G that contains u\V2 but not U2V\ weighs at least w{C). To prove this imagine contracting each 
edge x\X2 of G; call the resulting vertex x. Every vertex now has degree 2 in M. We will compute 
the weight of M by examining its edges in the contracted graph. An edge x±X2 £ M becomes a 
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loop at x of weight 0. An edge xy £ E \ E~ with 212/1,2:22/2 £ M or 212/2,222/1 G M becomes 2 
copies of xy, both with nonnegative weight. The other edges of M form cycles in the contracted 
graph. Each cycle is a cycle in G and so has nonnegative weight. One of the cycles contains edge 
uv, so it weighs at least w(C). Hence w(M) > w(C). 

Consider an even-cycle cover C that contributes to 5. ^ shows i 7^ j, i.e., u\Vi and U2^i occur 
with different multiplicities in C. The possibilities for {i,j} are {0, 1}, {0, 2}, and {1, 2}. C decom- 
poses into 2 perfect matchings. In all three cases one of the matchings of C contains exactly 1 of 
the edges U1V2, U2V1. That matching weighs at least w{C). The other matching has nonnegative 
weight, so C weighs at least w(C). In other words term*(<5) > w{C). □ 

Using this lemma with the scheme of Algorithm [7J gives the following. 

Algorithm 10 Checks whether a shortest cycle in undirected graph G with negative weights has 
weight < c. 

1: For = (1,2), (2,1), let d x P' 3 be the routine given by the Baur-Strassen theorem to 

compute the matrix of partial derivatives g— - — teim.y +nW (X/iLo 2/*) ( det(A(G)y w ) ) . 

2: Generate a random substitution a : X — > Z p for a prime p of order 0(n 4 ). 
3: For = (1, 2), (2, 1) compute the matrix <5 4J = d x f l ^\ a . 

4: Compute the matrix 6' with S' uv = x UliV2 5 u '% — x U2jVl Su,v 

5: Return true if 5' has a non-zero entry. 

Again a binary search gives 

Theorem 43. Let G = (V, E, w, W) be a weighted undirected graph without negative weight cycles. 
The weight of a shortest cycle can be computed in 0(Wn u ) time, with high probability. 



Algorithm 10 with c equal to the shortest cycle weight gives an edge uv with 8' uv 7^ 0, i.e., uv is 
on a shortest cycle. So a minimum weight perfect matching on G — u\, Vi corresponds to a shortest 
cycle. Hence we can state 

Theorem 44. Let G = (V, E, w, W) be a weighted undirected graph without negative weight cycles. 
A shortest cycle in G can be found in 0(Wn u ') time, with high probability. The same holds for a 
shortest st-path, for any given vertices s,t. 

10 Vertices Lying on Short Cycles 

For undirected graphs we only need to change the output of our algorithms: For Algorithm [7] (when 



there are no negative edges) or Algorithm 10 (in the general case) we find the set of vertices lying 



on cycles of length < c by changing the last step so that it returns {v S V : 3„ ue e S' vy _^= 0} . 
For directed graphs we apply the Y^a^q V % multiplication technique to Algorithm ^1 
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Algorithm 11 Computes the set of vertices lying on cycles of length < c in directed graphs. 
1: Let d x f be the routine given by the Baur-Strassen theorem to compute the matrix of partial 



derivatives 



a 



term 



c+nW 



dx u ,v VVJi "*?/ 
Generate a random substitution a : X 
Compute the matrix 5 = d x f\ a . 
Return {v G V : %,u)&E $(v,u) + 0}. 



,nW 



det((A0) + I)y 



Z p for a prime p of order 6(n ). 



Thus we obtain: 

Theorem 45. Let G be a weighted directed or undirected graph with integral weights in [— W, W] 
and no negative cycle. For any c the set of vertices lying on cycles of length < c can be computed 
in 0(Wn w ) time, with high probability. 
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